Fidelity in topological quantum phases of matter 
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Quantum phase transitions that take place between two distinct topological phases remain an unexplored area 
for the applicability of the fidelity approach. Here, we apply this method to spin systems in two and three 
dimensions and show that the fidelity susceptibility can be used to determine the boundary between different 
topological phases particular to these models, while at the same time offering information about the critical 
exponent of the correlation length. The success of this approach relies on its independence on local order 
parameters or breaking symmetry mechanisms, with which non-topological phases are usually characterized. 
We also consider a topological insulator/superconducting phase transition in three dimensions and point out the 
relevant features of fidelity susceptibility at the boundary between these phases. 
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I. INTRODUCTION 

In recent years topological phases have been intensively 
studied in condensed matter systems. These exotic states ap- 
pear in different contexts, such as fractional quantum Hall 
physics, spin liquids and topological insulators tH 0. HI B HI 
la, 01 • Their understanding is relevant to topological quantum 
computation, which provides the paradigm to store and ma- 
nipulate information in topologically non-trivial systems t^. 
In a seminal work ^ Kitaev introduced a spin model which 
can be exactly solved using a mapping to Majorana fermions 
coupled to a static Z2 gauge field ( |9], see also Ref. 1 10] for 
a perturbative approach). Generalizations of the Kitaev model 
with respect to lattice geometry, spatial dimension and local 
Hilbert space dimension have appeared recently [111, [T2I1 . 
In some regions of parameter space these models have non- 
trivial topological properties which cannot be described by 
any local order parameter. This prevents the direct applica- 
bility of the Landau-Ginzburg paradigm for the study of their 
phase transitions. Alternative ways of understanding critical- 
ity in systems with no local order parameter have been sug- 
gested in connection to quantum information. In particular, 
entanglement entropy (for a review see Ref. 1 13]) and the fi- 
delity approach {14, Oil [M \vjl fisll have attracted a lot of 



attention (for a review see Ref.|19l]). The reason why these 
two quantities can describe topological phases or determine 
their boundaries is due to the fact that they depend only on the 
properties of the ground state of the system and do not require 
a priori knowledge of any order parameter. 

In this paper we shall focus on the fidelity approach to topo- 
logical phase transitions ll20l I21I I22I1 . Much work has been 
done in understanding the nature of the phases in the Kitaev 
honeycomb model and, in particular, its fidelity has been stud- 
ied for the first time in Refs. [|23,.,24.1 . The present work aims 
to provide a more general understanding of fidelity in topo- 



logical phase transitions. For this porpuse we consider two- 
and three-dimensional extensions of the Kitaev honeycomb 
model. The models we chose differ from Kitaev's model in 
the geometry of the lattice, the interactions or the dimension 
of the ambient space and of the local Hilbert space. 



n. MODELS AND METHOD 

The dramatic changes in a many-body ground state | ^I^o ( A)) 
across a quantum phase transition can be captured by the fi- 
delity F = I (^'o(A)|^'o(A + SX)) I between two ground states 
corresponding to slightly different values of the set of param- 
eters defining the Hamiltonian H{{X}). Alternatively, one 
can use the fideUty susceptibility x = limsA^o ^^j,"^ ll25ll . 
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The fidelity quantifies how different two quantum states are. 
Given the drastic changes that take place across a quantum 
phase transition, one should expect the fidelity for nearby 
states to exhibit a drop across the boundary phase, signaling 
an enhanced distinguishability. Indeed, while in the thermo- 
dynamic limit the fidelity (susceptibility) has a drop (diver- 
gence) at the quantum critical point, for finite size systems 
this behavior is translated in the fidelity susceptibility being 
extensive away from criticality, while superextensive at crit- 
icality when the operator driving the transition is sufficiently 
relevant ifTsll . Moreover, a scaling analysis of the fidelity sus- 
ceptibility can be performed to extract the critical exponent ly 
of the correlation length ^. 

Here we apply this tool from quantum information theory 
to study the quantum phase transition between two distinct 
topological phases in three models: mosaic models, Kitaev's 
honeycomb model with an external magnetic field, and a 3D 
model exhibiting a transition between two distinct non-abelian 
topological phases. Let us first review them very briefly. 

An elegant way to look at the Kitaev model as a special 
instance of a more general class of models is provided by the 
mosaic classification lfT2ll . Mosaic models are defined on two- 
dimensional Bravais lattices constructed with a trivalent ver- 
tex building block, constrained by translational and rotational 
symmetry. The trivalent vertex is the local border between 
three polygons (see FiglT]). The number of edges of the three 
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(3-12-12) model (4-8-8) model 



FIG. 1: The trivalent vertex and the three different mosaic models 
constructed with it, which are considered in the present work. 




FIG. 2: The diamond lattice used in the three-dimensional extension 
of the Kitaev model. 

polygons (ei, 62, 63) specifies the mosaic model and only four 
possible such lattices are allowed: (6 — 6 — 6), corresponding 
to the Kitaev model; (3—12 — 12), corresponding to the model 
studied in lii, (4 - 6 - 12) and (4 - 8 - 8). 

The Hamiltonian of all mosaic spin models is given by 

u=x,y,z (i,j)eS(u) 

where S{u) is the set of edges in the w-direction. At every site 
the Hilbert space associated with the spin- 1/2 particle is iden- 
tified with a two-dimensional physical subspace of a four- 
dimensional Fock space M., introduced to solve the model. 
The Pauli operators acting on M. are represented by four Ma- 
jorana operators on M: — ib^c, — ib^c and — ib^c, 
with = 1, a/3 = -/3a, for a, /3 € {b'' ,by ,b^ ,c}, and 
a 7^ /3. The physical subspace is obtained through a projec- 
tion: 1^) e M D\i}) = with D = ib^'byb^c. The 



Hamiltonian ([T]l can be rewritten as = j 

with the operators Aij = iJuZij, where Zij = and 

(i, j) € S{u) {Aij — otherwise). The operators Zij com- 
mute with each other and with the Hamiltonian, so the space 
of the multi-spin system C can be decomposed as a direct 
sum C — 02 ^Zi where each sector is indexed by a set of 
eigenvalues {zy = ±l|(i,j) g S{u),u = x, ?;,z}ofthe 
operator Zij. Within each sector, the Hamiltonian reduces 
to a quadratic form in the Majorana operators a. The cou- 
plings of the Hamiltonian depend on a choice of vortex con- 
figuration, given by the eigenvalues of the plaquette operators 
Wp(„,) = - riijGapCn.) ' ^^^^ °f isogons 

{dp{ni) is the set of links belonging to plaquette p{ni )). 
The operators Wp(„.) commute with each other and with the 
Hamiltonian and since M^p^^ j = 1, their eigenvalues are 
Wp = ±1. A plaquette with Wp = +1 is a vortex free plaque- 
tte, while Wp — —1 corresponds to a vortex. In the following, 
we restrict ourselves to the vortex free sector, with a choice 
of Zij that preserves translational symmetry. T he g round state 
energy is at a minimum for this configuration f2M . 

Another modification of the Kitaev model that preserves its 
topological nature is the addition of a three-spin interaction of 
the form 

H = - E ^^j-^E^'^K' (2) 

u=x,y,z (jj')GS(u) i,j,k 

where the K-term is obtained through a perturbative expan- 
sion of a weak (Zeeman) magnetic field V = h ■ a iIztIi . In 
this case K ^ " jl ' and this Hamiltonian is assumed to ap- 
proximate the one with a Zeeman term when K <^ J, with 
Ju = J- The inclusion of this magnetic perturbation to the 
original model generates a topological phase with nonabelian 
anyonic excitations in the regime \Jz\ < | J^; | + | Jj, | , \Jy\ < 
\Jz\ + \ Jx\, \Jx\ < I 1 + I 1 ; while outside this region the 
excitations remain abelian. 

Recently analogues of the Kitaev model on three- 
dimensional lattices have been constructed ifTll |28|] . In the 
present work we focus on the model introduced in Ref. ifTlll 
since it presents quantum phase transitions between distinct 
non-trivial topological phases that, to our knowledge, have 
never been investigated with the fidelity approach. The model 
is defined on a diamond lattice and has a four-dimensional 
local Hilbert space (see Fig|2]). The Hamiltonian is the fol- 
lowing 

3 

^^-E-^" E +^^1)' (3) 

u=o (»j)es(u) 

where a and r are Pauli spin matrices associated to the two lo- 
cal spin- 1/2 degree of freedoms at each point in the lattice. As 
shown in ifTTIl Hamiltonian ([3]) can be mapped to a free Ma- 
jorana Hamiltonian which, in its ground-state sector, is given 
by 

ff = ^^J„ Y (A|At + AfA|), (4) 

«=0 (jj)eS(n) 
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and A^'^ are Majorana fermions introduced in the represen- 
tation of the PauH spin matrices (see Ref. Ii 1 lil for details), 
analogously to the previous model Hamiltonians. The addi- 
tion of next-nearest-neighbour interactions into (|4| allows for 
the presence of distinct non-trivial topological phases in the 
phase diagram, as we shall see in the next section. 

The spin realization of all the above models can be mapped 
onto a free Majorana fermion Hamiltonian 



(5) 



where Aij encodes the lattice structure and the parameters 
of the model. Once a unit cell of size s has been chosen the 
Hamiltonian can be written in direct space as 



HiA) 



E 

n,z-';m,/j. 



(6) 



where n and m are the positions of the unit cells in the lattice 
and /i and i' are the positions of the vertices inside the unit 
cell. Due to translational invariance A depends only on fi 
and m — n. Using a Fourier transform, we obtain 



(7) 



e *'"cr,a a complex fermion, N = 



withi'„(k) = i/v]vx;. 

sL** the system size, L the number of unit cells in one direc 
tion of the d-dimensional lattice and s the size of the unit cell, 
with a e {1, 2, .., s}. Then in general we have to deal with 
a free-fermion model characterized by s different bands. The 
ground-state is obtained filling the Fermi sea, where all the 
levels with negative energy e;3(k) < are occupied. 



ivi/o)^n^^(k)io) 



(8) 



/3,k 



if e;3(k) < 0. blili) = *t(k) . V/3(k), where V/3(k) is 
the eigenvector associated to the e;3(k) eigenvalue of 7i(k). 
The matrix J7(k) diagonalizing 7Y(k) has as column vectors 
V/5(k). The spectrum of the single-body Hamiltonian 7i(k) is 
symmetric around zero. 

The k-component of the many -body groundstate |5'o) = 
Ilk l^o)k in first quantization is given by a functional Slater 
determinant. The j-th component of the eigenstate corre- 
sponding to the /9-th band is represented by Vj^p (k) . Denoting 
with Sn the number of negative single particle energy bands, 
the first quantized wave function at fixed k is the Slater deter- 
minant of Sn particles that can occupy s„ bands 



(jl,---,Js„l*o)k = 



(9) 

withji G {1, 2, s}. Forthese particular models, the fidelity 
corresponds to the product over all k of the absolute value of 



the overlap of two Slater determinants at different parameter 
values ||(*o|*o)|| = 0k II (^o|*o)kll with 



(*o|*o>k=-^ E 



^.■=„4(k) ••• V,-.„,.„(k) 
VnA^) ■■■ ^.i,.„(k) 

^.■=„,i(k) ••• ^,.„,.„(k) 



.(10) 



The previous expression, using the properties of determi- 
nants, can be rewritten as (for a proof see pag.29 1 of Ref. i29ll ) 



Vi(k).Vi(k) ••• Vi(k).V,„(k) 



V, (k).Vi(k) ••• V, fk).V,„(k) 



(11) 



The above formula will be used in the evaluation of fidelity 
for the models that we consider. 



III. RESULTS 

In this section we study how the fidelity susceptibility be- 
haves for the mosaic models (4 — 8 — 8) and (3 — 12 — 12), 
for the model dU, and for the three-dimensional Kitaev model 
(0. 



A. Mosaic models 

Let us first focus on the (4—8 — 8) mosaic model ||3C 
Choosing a four-site unit cell, a Fourier transform of the 
Hamiltonian gives 



i7=i^*t(k)H(k)*(k), 



(12) 



where 



w(k) = ( r"., t V (13) 

— (°^k 1' '^k 2' ^k 3' °^k 4) Fourier trans- 

form of the Majorana operators {'^i}f^i, ct — 
diag[exp (— ifc2), — exp (z/ci)], ki ~ k ■ ni, A:2 = k • n2, 
ni = (1,0) and n2 = (0, 1). The system size is TV = 2L2. 
This model presents a quantum phase transition between two 
gapped phases with abelian anyons when J| = + Jy. The 
phases are algebraically distinct, though they can be related 
by rotational symmetry. 

For the (3—12 — 12) mosaic model ll6l [3lll . a Fourier trans- 
form with a six-unit cell renders the Hamiltonian in the form 
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I 4. 4. 



— ikni t 

"k,l 



ak,2 



''k,4 



Ok, 5 



''k,4'' 



ak,6 - aicsftk^e) + he, 



(14) 



with m = (1/2, V3/2) and na = (1/2, -V3/2). The single 
particle Hamiltonian A(k) is 6 x 6, and the system size is 
N = 6L^. This system has a quantum phase transition at 
J — \/3 from a topological phase with abelian anyons to one 
with nonabelian anyons, spontaneously breaking time reversal 
symmetry, without the need of an external magnetic field. 

The numerical results for the fidelity susceptibility x = 
lim^A-^o — 21nF/(5A^ for the (4 — 8 — 8) mosaic model are 
shown in Fig. |3]for different system sizes, taking — Jy ~ 
l,Jz G (1.21, 1.61). We see that x is an extensive quantity 
off criticality, while it is superextensive at the critical point 
= ^/2. 
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FIG. 4: (Color online), (a) Scaling analysis of the fidelity suscepti- 
bility at criticality for the (4 — 8 — 8) mosaic model, for system sizes 
L £ [300, 1000] in steps of 100, with 5 J = 10"''. We obtain that 
X /N ^ L", with fi = 0.1523±0.0001. (b) Data collapse for the 
fidelity susceptibility around the critical point — -\/2. All curves 
collapse when plotted as a function of the dimensionless quantity 
L^^"{J^ - =''')/ Jf, with critical exponent ly = 1.105 ±0.070 



To quantify the extent of the collapse, we use a method 
similar to the one described in Ref. 113311 . obtaining the value 
1/ = 1.105 ± 0.07. From there, we conclude that a = ^ ^ 
0.138 ± 0.009. The result for the data collapse is shown in 
Fig.|4](b). 

A similar analysis has been performed for the (3 — 12—12) 
mosaic model [6], confirming that the fidelity is superexten- 
sive at criticality, with /i ~ 0.155 ± 0.009 as seen on Fig. 

El 

However, we have not been able to obtain a satisfactory data 
collapse for the (3 — 12 — 12) mosaic model, probably due to 
numerical errors. 



FIG. 3: (Color online). Fidelity susceptibility x per site for the 
(4—8—8) mosaic model around the quantum critical point — \/2, 
separating two topological phases (Jx = Jy — 1). x shows exten- 
sive behavior off criticality, while it is superextensive at the critical 
point. 

We perform finite-size scaling for the susceptibility around 
the critical point = \/2- For a finite sample size we denote 
the point at which x is maximum as x™^^, evaluated at = 
J™'^'^. It scales Uke x""'''/N cx L^', with /i > 0, while in 
the thermodynamic limit x/iV c>c 1/| J' — Jz\"- In Fig- 111(a) 
we plot logx"^^^ vs log!/, for system sizes L £ [300, 1000] 
in steps of 100, with SJ = 10~^. We obtain superextensive 
scaling given by ^ = 0.1523 ± 0.0001. 

Furthermore, we perform a data collapse, using the scaling 
ansatz lHai 



xiJ) 

T^max 



/( 



T — 
V2 



(15) 



B. Kitaev model with 3-spin interaction 

This model can be exactly diagonalized using the tech- 
niques we have summarized above. In particular, the Hamil- 
tonian in equation ^ takes the form H = \ ^ ■ ^ AijCiCj, 

with Aij — iJuZij + KJ2k ^ikZjk, and for the vortex-free 
case, the spectral matrix in momentum space is 



5(k) i/(k) \ 

pikni 



(16) 



where /(k) = Jz + J^e"^"^^ + J2e*'^"=, .g(k) = 
2if(sink • rii — sink • n2 + sink • (n2 — ni)), with rii = 
(1, 0) and n2 = (0, 1). The system size is iV = 21"^. 

The inclusion of an external magnetic field gives rise to a 
transition between phases with abelian and nonabelian anyon 



5 



ln(x )=2.1551n(L)-0.928 



6.2 6.4 

ln(L) 



^ 15 
14.5 
14 



-(a) 


(b) 


/ ^ 


- 






/ \ 




itiyX J— z. 1 / jjirnLj+u.yzUo 
of 


/ 


=° * 

8° 

* % 


- 




8 




s 

I 
% 








Ip 




f 




&^ 


r 1 . 1 < 1 < 1 < 1 < 1 . 




.1,1 


, 1 4- 



5.8 6 6.2 6.4 6.6 6.8 7 -0.6 -0.3 0.0 0.3 0.6 

ln(L) -J "^'')/J ' 



?5 



FIG. 5: (Color online). Scaling analysis of the susceptibility at crit- 
icality for the (3 - 12 - 12) mosaic model, for L € [301, 1001] in 
steps of 100. The fidelity signals the quantum phase transition sep- 
arating the two distinct topological phases at J = \/3, behaving as 
^maxy^ - L", with fi = 0.155 ± 0.009 



excitations, contained in the region | | < | J^^ | + | Jj, | , | Jj, | < 



J. 



JJ, JJ < J„ 



\Jz \ of the parameter space. 



The numerical results for yf^^^ /N cx for the honey- 
comb model with external magnetic field are shown in Fig. 
|6](a), where we obtain ^ = 0.1735 ± 0.0001, signaling su- 
perextensive scaling at criticality. The corresponding data col- 
lapse is shown in Fig. |6](b), with v = 1.10 ± 0.05 for the 
critical exponent of the coiTelation length, and a ^ [ijv ^ 
0.158 ±0.007. 

For the models we have considered so far we see that the 
fidelity susceptibility is able to detect the phase transition be- 
tween distinct topological phases, besides providing useful in- 
formation about the critical point itself (i.e. the exponent v of 
the correlation length). 



FIG. 6: (Color online), (a) Scaling analysis of the susceptibility for 
the Kitaev model with a weak Zeeman magnetic field, at the critical 
point = 1 /2. The system size range is L G [301 , 1001] in steps of 
100, and the coupling K = 1/15. The data implies that x™'"'/A'' oc 
i'', with ii = 0.1735 ± 0.0001 at criticality. (b) Data collapse of 
the fidelity susceptibility x around the critical point = 1/2 ( Ja; = 
3y = 0.5(1 - J^)) and K = 1/15. We obtain v = 1.10 ± 0.05 for 
the critical exponent of the correlation length. 



removes the degeneracy ifTlll . In momentum space the Hamil- 
tonian ^ with the above perturbations can be written as 

H = iY,{a\, at,,bi„bt,)'H{k){alalbl hi), (19) 

k 

with Ti.{k) given by 



, e(k) ^m 
^^^> ^ [ -z$*(k) -e(k) 



(20) 



Defining the three component vectors sq = l/4(— 1, 1, —1), 
si = 1/4(1,1,1), S2 = l/4(-l,-l,l) and S3 = 
1/4(1, —1, —1) (see Fig|2]l, the blocks of the matrix (l20l l are 
given by 



C. 3 dimensional model 

The Hamiltonian in equation ^ has lines of zeros in 
momentum space with a vanishing gap. The addition of 
next-nearest-neighbour interactions that preserve time rever- 
sal symmetry of the form 



jk 



(17) 



u=0 



e(k) = e^(k)cr^ + e"(k)CT^ 



e^(k) =i^sin^^^-^ 



(21) 



(22) 



Sill + Sin — h sm 



(23) 
(24) 



jk 



(18) with (t'' the two-dimensional identity matrix and , z the 
usual Pauli matrices. 
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When all couplings J are equal, there are still three de- 
generate points in the Brillouin zone at (27r, 0, 0), (0, 27r, 0) 
and (0, 0, 2tt). These degeneracies can be made massive by 
adding a small anisotropy in one of the couplings J. Follow- 
ing [llj we set — 1, ^2,3,4 — 2 and Ji = 2 + 6Ji. The 
free parameters of the model are then and 5Ji. The phase 
diagram presents three topologically distinct phases, distin- 
guished by the winding number For positive 5Ji and positive 
tz the winding number is +1; for positive 5Ji and negative tz 
the winding number is —1, and it is zero in the negative 6Ji 
region 1 1 1]. There are two topological phase transitions at the 
two critical lines 6Ji = and = 0. 

We have studied this model using the fidelity approach and 
in the following we present the results for the scaling of the 
fidelity susceptibility around the two phase transitions of the 
model. In Fig|7](a) we show the log-log plot of the fidelity sus- 
ceptibility at the tz = critical line, while Fig|2](b) presents 
analogous results for the transition across the 5Ji ~ Q line. 
In both cases the fidelity susceptibility has a peak at the crit- 
ical point. The scaling analysis suggests that the fidelity sus- 
ceptibility is extensive within the numerical precision that we 
could obtain from a first set of data. However, the numerical 
accuracy of this result cannot exclude a weak super-extensive 
scaling at the critical point. A more accurate analysis of this 
model will be presented in a future work. 
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FIG. 7: (a) Log-log plot of the fidelity susceptibility at the critical 
line for 5Ji = 1, and linear system size L G [200, 300] in steps of 
20. (b) Log-log plot of the fidelity susceptibility at the critical line 
for tz = 0.1, and linear system size L G [200, 300] in steps of 20 



IV. CONCLUSIONS 



The fidelity approach proves to be a useful tool to detect 
and characterize quantum phase transitions which take place 
between distinct topological phases. We have applied it to 
the family of mosaic models which presents a transition be- 
tween topological phases with abelian and nonabelian anyons, 
as well as the Kitaev model with an external magnetic field, 
where a transition from an abelian to a nonabelian phase takes 
place. We were able to show superextensive scaling of the 
fidelity susceptibility at criticality, characterizing its behavior 
in the thermodynamic limit. We also extracted the critical ex- 
ponent of the correlation length for these models. Moreover, 
we studied the transition between a topological insulator and 
a topological superconductor in a three-dimensional diamond 
lattice, where the method still detects the boundary between 
different phases. In this case, superextensivity of the fidelity 
susceptibility turns out to be more difficult to be established. 
This is probabily due to the increased dimensionality which 
weakens the features of the quantum phase transitions. 

We have thus extended the applicability of the fidelity ap- 
proach to the boundary that separates different topological 
phases and proved its usefulness in this context as well. Since 
it is not always clear how to detect and distinguish differ- 
ent topological phases of matter we believe it is important to 
check the validity of a method based solely on the geometri- 
cal properties of the groundstate, and that does not require the 
introduction of an order parameter. 
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